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The Sagdeev potential technique has been considered to find the necessary condition 
for the existence as well as the polarity of dust ion acoustic solitary structures at 
M = M c , in an unmagnetized nonthermal plasma consisting of negatively charged 
dust grains, adiabatic positive ions and nonthermal electrons, where M c is the small- 
est possible value of the Mach number M. Depending on the analytical theory, a 
compositional parameter space showing the nature of existence of solitary structures 
at M = M c has been presented. The compositional parameter space clearly indicates 
the upper bounds of the nonthermal parameter (/3j) associated with the nonthermal 
distribution of electrons for the existence of negative and positive potential solitons 
at M — M c . The compositional parameter space also shows that the present system 
does not support the coexistence of both negative and positive solitons at M — M c . 
A critical value of the nonthermal parameter (Pi = fi c ), which acts as sink of negative 
solitons and the source of positive solitons at M — M c , has been found. However, 
no solitary wave is possible at (3\ = (3 C . The compositional parameter space again 
shows that the present system supports the negative potential double layer along the 
curve Pi = [3i c in /^/Si-parametric plane provided that /i > fi r , i.e., for each \i lying 
within the interval /j, r < \l < 1, there exists one and only one Pi c of pi such that one 
can get negative potential double layer at M = M c when Pi = /3 lc , where fi is the 
ratio of unperturbed number density of nonthermal electrons to that of ions. From 
the compositional parameter space, it has been observed that the negative potential 
solitons at M = M c are separated by the curve Pi = Pi c . A finite jump between 
amplitudes of negative potential solitons separated by the curve pi = Pi c has been 
observed. The effect of the parameters on the amplitude of the solitary structures 
have been discussed. 
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I. INTRODUCTION 



In the last few decades dusty plasma attracts the researchers most, not only because of 
different types of wave propagation in dusty plasma due to their involvement in the study 
of astrophysical and space environments but also for its great variety of new phenomena as- 
sociated with waves and instabilities.- - - Dusty plasmas essentially consists of fairly massive, 
charged particulates in addition to electrons and ions in the usual two component plasma. 
Acoustic wave modes in dusty plasma have received a great deal of attention since the last 
decade. Depending on different time scales, there can exists two or more acoustic waves in 
a typical dusty plasma. Dust Acoustic (DA) and Dust Ion-Acoustic (DIA) waves are two 
such acoustic waves in a plasma containing electrons, ions, and charged dust grains. 

Shukla and Silin^ were the first to show that due to the quasi neutrality condition n e o + 
ridoZd = n i0 and the strong inequality n e0 <C n iQ (n e0 , n i0 , and na are, respectively, the 
number density of electrons, ions, and dust particle, where Zd is the number of electrons 
residing on the dust grain surface), a dusty plasma (with negatively charged static dust 
grains) supports low-frequency DIA waves with phase velocity much smaller (larger) than 
electron (ion) thermal velocity. In case of long wavelength limit the dispersion relation of 
DIA wave is similar to that of IA wave for a plasma with n e o = n^o and Tj <C T e , where Tj(T e ) 
is the average ion (electron) temperature. Due to the usual dusty plasma approximations 
(n e0 <C n i0 and T, t ~ T e ), a dusty plasma cannot support the usual Ion-Acoustic (IA) waves, 
but the DIA waves of Shukla and Silin^ can. Thus DIA waves are basically IA waves, 
modified by the presence of heavy dust particulates. The theoretical prediction of Shukla 
and Silin^ was supported by a number of laboratory experiments.— - — The linear properties 
of DIA waves in dusty plasma are now well understood.— 1 ^ - — 

Due to a balance between the wave dispersion and nonlinearity, there exists the possibility 
of DIA solitary waves (DIASWs) in a dusty plasma. The DIASWs studied using the reduc- 
tive perturbation technique which gives Korteweg-de Vries (KdV) equation or Korteweg- 
de Vries Burgers (KdVB) equation as nonlinear evolution equations.— Bharuthram and 
Shukla^ studied the nonlinear waves associated with DIA waves. Mamun and Shukla^ 
have investigated the condition for existence of positive and negative potential DIASWs, 
and examined how the ion-fluid temperature modifies the amplitude of DIASWs. In pres- 
ence of negatively charged dust grains, one can have DIASWs with negative potentials, when 
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the electron number density is less than one third of the positive ion number density^ 1 ^. 
Verheest et aP^ have shown that in the dust-modified ion acoustic regime, negative struc- 
tures can also be generated, beside positive potential soliton for electrons with a polytropic 
index 7 e ^ 1, as well as for Boltzmann electrons. The effect of ion-fluid temperature on 
DIASWs structures have been investigated by Sayed and Mamurt^ in a dusty plasma con- 
taining adiabatic ion-fluid, Boltzmann electrons, and static dust particles. In the above 
mentioned works the nonlinear properties of DIASWs investigated for M > M c . 

In most of the earlier works, Maxwellian velocity distribution function for lighter species 
of particles has been used to study DIASWs and DIA double layers (DIADLs). However 
the dusty plasma with nonthermally/suprathermally distributed electrons observed in a 
number of heliospheric environmenta^^Sr— . Therefore, it is of considerable importance to 
study nonlinear wave structures in a dusty plasma in which lighter species (electrons) are 
nonthermally/suprathermally distributed. Recently Baluku et alr^ investigated DIASWs 
in an unmagnetized dusty plasma consisting of cold dust particles and kappa distributed 
electrons using both small and arbitrary amplitude techniques. 

In the present investigation we have considered the problem of existence as well as the po- 
larity of DIASWs and DIADLs in a nonthermal dusty plasma consisting of negatively charged 
dust grains, adiabatic positive ions and nonthermal electrons at the smallest possible value of 
the Mach number M. The Sagdeev potential approach^ has been considered to investigate 
the existence of solitary wave and double layer at M — M c , where M c is the lower bound of 
the Mach number M, i.e., solitary wave and/or double layer solutions of the energy integral 
start to exist for M > M c . In most of the earlier works, solitary wave and/or double layer so- 
lutions have been investigated for M > M c . However, some recent investigations^ 1 ^— have 
shown that finite amplitude solitary wave can exist at M = M c in the parameter regime 
where solitons of both polarities exist. The numerical observations^^— of the solitary 
wave solution of the energy integral at M = M c , motivate us to investigate the existence of 
DIASW and DIADL. However, we have approach in a different way without considering the 
parameter regime of the system for M > M c . Three basic parameters of the present dusty 
plasma system are /i, a and px, which are respectively the ratio of unperturbed number 
density of nonthermal electrons to that of ions, the ratio of average temperature of ions 
to that of nonthermal electrons, a parameter associated with the nonthermal distribution 
of electrons. Nonthermal distribution of electrons becomes isothermal one if (3x = and 
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consequently, for isothermal electron species, the present dusty plasma contains only two 
basic parameters /i and a. It has been found that solitary wave or double layer do not 
exist for < \i < \i v . For isothermal distributed electrons, PPSWs exist for \i p < // < /i c , 
whereas for all /i lies within fi c < \x < fir, there exist NPSWs. At \i — fi c , no solitary 
wave or double layer solution is possible. However, there do not exists any double layer 
solution for isothermal electrons. For nonthermal distributed electrons, we have found two 
critical values of (3\, viz., (3\ a and (3 C are the upper bounds of (3\ for the existence of PPSWs 
and NPSWs, respectively. However, the coexistence of both NPSWs and PPSWs are not 
possible in the present system. Again, at f3\ = (3 C , no solitary wave or double layer solution 
is possible. We have also found the existence of negative potential double layer (NPDL) in a 
portion of the solution space. Depending on the nature of solutions of DIA waves, we have 
fi p , /i c , and fi r such that the entire parametric space can be delimited into three sub intervals 
of /I, viz., /i p < ji < fi c , fi c < ji < fi r , and fi r < \i < ht- For each // in \i r < \i < fix, 
there exists a value j3\ c of fix, such that one can found the existence of NPDL. There is 
an interesting phenomena occur around the point /?i = /3i c . There exists a finite jump in 
amplitude of NPSWs at 0x = Pic, the jump being the amplitude of the double layer exists at 
0x — Pic- The amplitude of NPSW decreases with increasing f3\ and ultimately, demolished 
at 0x = Pc- On the other hand, PPSW creates from /?i = /3 C and the amplitude of PPSW 
increases with f3\, having maximum amplitude at /?i = ftxa- Throughout the investigation 
we have considered [it as the upper bound of fi, which is strictly less than 1, since for jj, — > 1, 
the effect of charged dust on DIA wave is negligible. 

The present paper is organized as follows: In Sec. [Til the basic equations are given. 
The energy integral and the Sagdeev potential has been constructed in Sec. IIHI Physical 
interpretation of the energy integral for the existence of solitary wave and double layer 
solutions is also given in this section. In Sec. IIVI an analytical approach has been presented 
to determine the polarity of the solitary waves. The analytical approach has been verified 
numerically in Sec. |Vj Finally, a brief conclusion has been given in Sec. IVIl 

II. BASIC EQUATIONS 

The governing equations describing the nonlinear behavior of DIA waves, propagating 
along x-axis, in collisionless, unmagnetized dusty plasma consisting of negatively charged 
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immobile dust grains, adiabatic positive ions and nonthermal electrons are the following: 

Ui,t + (riiUi) x = 0, (1) 



u i)t + UiU i)X = -<f) x - an i 1 p itX , 



(2) 



Pi,t + Uip ijX + "YPiU iy 



i,x 



(3) 



XX 



(n e0 /n i0 )n e -m + (Z d n d )/n i0 , 



(4) 



where the parameter a = Ti/T e . 

Here we have used the notation q or (ipi) q for dipi/dq and iii, n e , Ui, Pi, <f), x and t are, 
respectively, the ion number density, electron number density, ion velocity, ion pressure, elec- 
trostatic potential, spatial variable and time, and they have been normalized by (unper- 
turbed ion number density), n e o (unperturbed electron number density), q(= a/ (K B T e )/mi) 
(ion-acoustic speed), n i0 K B Ti, K B T e /e, A jDem (= a/ (K B T e )/ (47m i0 e 2 ))(Debye length), and 
oo^ i 1 (= >/mi/ ' (Airnioe 2 )) (ion plasma period). Here 7(= 3) is the adiabatic index, K B is the 
Boltzmann constant, Tj and T e are, respectively, the average temperatures of ions and elec- 
trons, rrii is the mass of an ion, n d is the dust number density, Z d is the number of negative 
unit charges residing on dust grain surface, and e is the charge of an electron. 

The above equations are supplemented by nonthermally distributed electrons as pre- 
scribed by Cairns et alM- for the electron species. Actually, in a number of heliospheric 
environments, dusty plasma contains nonthermally distributed ions or electrons. Therefore, 
it is of considerable importance to study DIASWs and DIADLs in dusty plasma in which 
electrons are nonthermally distributed. Nonthermal distribution of any lighter species of 
particles (as prescribed by Cairns et alM- for the electron species) can be regarded as popu- 
lation of Boltzmann distributed particles together with a population of energetic particles. 
This can also be regarded as a modified Boltzmannian distribution, which has the property 
that the number of particles in phase space in the neighbourhood of the point v = is much 
smaller than the number of particles in phase space in the neighbourhood of the point v — 
for the case of Boltzmann distribution, where v is the velocity of the particle in phase space. 
This type of velocity distribution is often termed as Cairns distribution and was considered 
by many authors in various studies of different collective processes in plasmas and dusty 
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plasmas^"—. Under the above mentioned normalization of the dependent and independent 

variables, the normalized number density of nonthermal Cairns distributed electrons can be 
written as 

n e = (1 -/?!</> + (5) 

where 

Pi = (6) 

H l + 3ai' v ; 

with oi\ > 0. Here Pi and consequently a\ are the parameters associated with nonthermal 
distribution of electrons. From ([6]) and inequality ot\ > 0, it can be easily checked that the 
nonthermal parameter Pi is restricted by the following inequality: < Pi < 4/3. However 
we cannot take the whole region of Pi (0 < Pi < 4/3). Plotting the nonthermal velocity 
distribution of electrons against phase velocity {y ) of electrons in phase space, it can be 
easily shown that the number of electrons in phase space in the neighbourhood of the point 
v = decreases with increasing /3i and the number of electrons in phase space in the 
neighbourhood of the point v = is almost zero when Pi — > 4/3. Therefore, for increasing 
values of Pi distribution function develops wings, which become stronger as Pi increases, 
and at the same time the center density in phase space drops, the latter as a result of the 
normalization of the area under the integral. So, we should not take values of Pi > 4/7 since 
that stage might stretch the credibility of the Cairns model too far—. So here we consider 
the effective range of pi as follows: < X < iT , where (3 1T = 4/7 ps 0.571429. When 
Pi = 0, nonthermal distribution of electrons becomes isothermal one. 

Now introducing a new parameter /i = n e0 /nio, the charge neutrality condition, 

ZdJido + n e0 = n i0 , (7) 
can be written as 

= 1 - V, (8) 

so that the Poisson equation (T4J) can be written as 

4>xx = l^n e - rii + (1 - //). (9) 

We note from ([S]) that 1 — fi must be greater than zero, i.e., < /i < 1. When /i — > 1, the 
effect of negatively charged dust grains on DIA wave is negligible and so, we restrict \i by 
the inequality < \i < fir, where \it is strictly less than 1. 



III. ENERGY INTEGRAL AND THE SAGDEEV POTENTIAL 



To study the arbitrary amplitude time independent DIA solitary structures we make all 
the dependent variables depend only on a single variable £ = x — Mt, where the Mach 
number M is normalized by q. Thus in the steady state, ([I])-© and flU]) can be written as 

- Mn it s + (niUi)^ = 0, (10) 

- Mu it 5 + UiU it ^ = - fa - an^ l p it ? , (11) 

- Mp i: $ + UiPi, £ + 3piUi, 5 = 0, (12) 

<f>tf = fm e -ni + (l-(j). (13) 

Using the boundary conditions, 

rii -)■ I, pi — y 1, Ui ->■ 0, — y 0, 0^ ->■ as |£| ->■ oo, (14) 

and solving (ITDjl - ( TI2"]) . we get a quadratic equation for n 2 and the solution of this quadratic 
equation can be put in the following form: 
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n? 



(v / *WH - M - </>) > ( 15 ) 



where 

T/w , (M-^3a) 2 (M + v 7 ^) 2 

*(M) = ^ $(M) = ^ (16) 

From the expression of nf as given by (fTBT) we see that n 2 is nonnegative if and only if 



a/^(M) - and a/$(M) - both are real. Again y/ty(M) - and a/$(M) - both are 
real if and only if \?(M) — and $(M) — both are real and nonnegative. But \I/(M) — 
and $(M) — both are real and nonnegative if and only if < \&(M) and < $(M), 
respectively, and finally the conditions < ^f(M) and < $(M) hold simultaneously 
if and only if < Min{^(M), $(M)}. Now from ([11]) we see that *(Af) < $(M) (^> 
Min{\l/(M), $(M)} = \I/(M)) and consequently ffTBT) gives theoretically valid expression of 
n 2 if and only if < \I/(M). Therefore, for the present problem we are restricted by the 
inequality: < \1/(M) = (M — \/3a) 2 /2. But numerically it is not possible to get the 
correct value of from ( jl5j) for small value of a because the denominator of the right hand 
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side of (fl5|) is 6a. To remove a from the denominator of the right hand side of (fl5|) we 
first of all note that (a/^M) - + \f$(M) - 0) 2 > for all < *(M) and consequently 
we can multiply the numerator and the denominator of the right hand side of (fl5l) by 



(a/^(M) - + a/$(M) - 0) 2 and we get the following expression of rij. 

x/2M 



71; 



(17) 



From (TlTj) . we see that this equation gives both theoretically and also numerically correct 
expression of even when a = if < \1/(M). Now if we put a = in (flTl) . we get the 
same expression of rij for the cold dust fluid with $(M) = \I/(M) = 

Now integrating (fT3|) with respect to and using the boundary conditions (fT4|) . we get 
the following equation known as energy integral with V{4>) = V^(M, 0, /i, a) as Sagdeev 
potential or pseudo-potential. 

where 

7(M, 0, ft, /i, a) = V5 - /iK - (1 - /i)0, (19) 
K = (1 + 3A - 3A0 + /3x0 2 )e^ - (1 + 3^), (20) 

^ = M 2 + a-n i (M 2 + 3a-20-2an J 2 ). (21) 

Although V(4>) is a function of 0, M, and a, but for simplicity, we can omit any 

argument from V(M, 0, /3\, fi, a) when no particular emphasis is put upon it. For example, 
when we write V(<p) = V(M,(p), it is meant that the arguments fi and a assume fixed 
values in their physically admissible range and for the fixed values of /3i, fi and a, nature of 
solitary structures as obtained from the energy integral (fl8l) depends only on M. 

The energy integral ffl8|) can be regarded as the one-dimensional motion of a particle 
of unit mass whose position is at time £ with velocity d(p/dC, in a potential well ^(0). 
The first term of the energy integral can be regarded as the kinetic energy of a particle 
of unit mass at position and time £ whereas V(<ft) is the potential energy of the same 
particle at that instant. Since kinetic energy is always a non-negative quantity, V{<p) < 
for the entire motion, i.e., zero is the maximum value for V(<ft). Again from ffl8|) . we find 
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d 2 (f>/d^ 2 + V'(4>) = 0, i.e. the force acting on the particle at the position is —V'(4>). 
Suppose, ^(0) = ^'(0) = 0, therefore, the particle is in equilibrium at = because the 
velocity as well as the force acting on the particle at = are simultaneously equal to 
zero. Now if = can be made an unstable position of equilibrium, the energy integral can 
be interpreted as the motion of an oscillatory particle if V(4> m ) = for some m ^ 0, i.e., 
if the particle is slightly displaced from its unstable position of equilibrium then it moves 
away from its unstable position of equilibrium and it continues its motion until its velocity 
is equal to zero, i.e., until takes the value m . Now the force acting on the particle of 
unit mass at position = m is — V(0 m ). For m < 0, the force acting on the particle at 
the point = m is directed towards the point = if — V'((f> m ) > 0, i.e., if V'(4> m ) < 0. 
On the other hand, for m > 0, the force acting on the particle at the point = m 
is directed towards the point = if — V'(<p m ) < 0, i.e., if V'(4> m ) > 0. Therefore, if 
V(<f>m) > (for the positive potential side ) or if V"(0 m ) < (for the negative potential side 
) then the particle reflects back again to = 0. Again, if V((j) m ) = V'((p m ) = then the 
velocity d<j)/dl; as well as the force d 2 (f)/d£ s 2 both are simultaneously equal to zero at = m . 
Consequently, if the particle is slightly displaced from its unstable position of equilibrium 
(0 = 0) it moves away from = and it continues its motion until the velocity is equal 
to zero, i.e., until takes the value = m . However it cannot be reflected back again at 
= as the velocity and the force acting on the particle at = m vanish simultaneously. 
Actually, if V'((j) m ) > (for m > 0) or if V'((f) m ) < (for m < 0) the particle takes 
an infinite long time to move away from the unstable position of equilibrium. After that 
it continues its motion until takes the value m and again it takes an infinite long time 
to come back its unstable position of equilibrium. Therefore, for the existence of a PPSW 
(NPSW) solution of the energy integral f|T8|) . we must have the following: (a) = is the 
position of unstable equilibrium of the particle, (b) V(4> m ) = 0, V'(4> m ) > (V(0 m ) < 0) 
for some m > (0 m < 0), which is nothing but the condition for oscillation of the particle 
within the interval min{0, m } < < max{0, m } and (c) V(0) < for all < < m 
(0 m < < 0), which is the condition to define the energy integral ( |T8l) within the interval 
min{0, m } < < max{0, m }. For the existence of a PPDL (NPDL) solution of the energy 
integral (|T8|) . the conditions (a) and (c) remain unchanged but here (b) has been modified 
in such a way that the particle cannot be reflected again at = 0, i.e., the condition (b) 
assumes the following form: V((p m ) = V'((p m ) = 0, V"((f> m ) < for some m > (0 m < 0). 
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The above discussions for the existence of solitary waves and double layers are valid if 
= is an unstable position of equilibrium, i.e., if V"(0) < along with V(0) = V'(0) = 0. 
In other words, = can be made an unstable position of equilibrium if the potential energy 
of the said particle attains its maximum value at = 0. Now, the condition V"(0) < gives 
a lower bound M c of M, i.e., V"(0) < <S> M > M c , V"{0) > <£> M < M c , and 
V"(0) = O M = M c o V"(M C , 0) = 0. This M c is, in general, a function of the 
parameters involved in the system, or a constant. Therefore, if M < M c , the potential 
energy of the said particle attains its minimum value at = 0, and consequently, = is 
the position of stable equilibrium of the particle, and in this case, it is impossible to make 
any oscillation of the particle even when it is slightly displaced from its position of stable 
equilibrium, and consequently there is no question of existence of solitary waves or double 
layers for M < M c . In other words, for the position of unstable equilibrium of the particle 
at = 0, i.e., for V"(0) < 0(<^ M > M c ), the function V((f>) must be convex within a 
neighborhood of = and in this case both type of solitary waves (negative or positive 
potential) may exist if other conditions are fulfilled. Now suppose that V"(M C , 0) = and 
also V"'(M C , 0) = 0, then if V""(M C , 0) < 0, the potential energy of the said particle attains 
its maximum value at = and consequently, = is the position of unstable equilibrium. 
On the other hand if V"(M c ,0) = 0, V"'(M c ,0) = 0, and V""(M c ,0) > 0, the potential 
energy of the said particle attains its minimum value at = and consequently, = 
is the position of stable equilibrium of the particle and in this case there is no question of 
existence of solitary wave solution and/or double layer solution of the energy integral ( TT8l) . 

However if V(M C , 0) = V'(M e , 0) = V"(M e , 0) = and V"'{M C , 0) < 0, then with the 
help of simple mathematics, it can be shown that V(M C , 0) is a concave function in (— e , 0) 
and V(M C , 0) is a convex function in (O,0 e ), where e is strictly positive quantity but 
can be made sufficiently small. The out line of proof of this result is given in Appendix. 
Therefore, if a particle, placed at = 0, be slightly displaced towards the positive potential 
side it falls within the interval (0, <f> e ) and due to the convexity of the function V(M C , 0) 
within the interval (0, <f> e ), it moves away from = and it continues its motion until its 
velocity is equal to zero, i.e., until takes the value = m > 0, where V(M C , m ) = 
and in this case if V'(M C , m ) > (V (M c , m ) = and V"(M C , m ) < 0), one may get a 
PPSW (PPDL) as a solution of the energy integral (fl8|) . On the other hand, if the same 
particle is slightly displaced from = towards the negative potential side, it falls within 
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the interval (— e , 0) and due to the concavity of the function V(M C , 0) within the interval 
(— e ,O), it moves towards the point = and consequently, there does not exist any 
solitary wave or double layer solution in the negative potential side. By similar arguments, 
we can conclude that if V(M, 0) = V'(M, 0) = V"(M e ,0) = and V"'{M C ,0) > then 
one may obtain a NPSW (NPDL) as a solution of the energy integral (1181) . whereas there 
does not exist any PPSW or PPDL. At this point, it is important to note that we can not 
expect the coexistence of both negative and positive potential solitary waves at M — M c . 
Therefore, the necessary conditions for the existence of PPSW or PPDL at M — M c are 
V(M C , 0) = V'(M e , 0) = V"(M C , 0) = and V"'(M C , 0) < whereas the necessary conditions 
for the existence of NPSW or NPDL at M = M c are V(M C , 0) = V'(M C , 0) = V"(M C , 0) = 
and V"'(M c ,0) > 0. 



IV. POLARITY OF SOLITONS AT M = M c 

For the present problem, it can be easily checked that V(M, 0) = V'(M, 0) = for any 
value of M and also for any values of the parameters /i and a, whereas the condition 
V"(M, 0) < gives M > M c , where 

^^^r^fty < 22) 

From (|22|) . we see that for M c to be real and positive, we must have > and < fli < 1. 
As the effective range of fi 1 is < ft < 1T , where (3 1T = 4/7 w 0.571429, M c is well-defined 
as a real positive quantity for all < /i < //^ and < 0i < ftr- 

From the discussions of the previous section we have seen that to determine the polarity 
of solitary wave or double layer at M = M c , the function V'"(M C , 0) plays an important role 
and for the present problem that function is given by 

V"'(M C , 0) = fi{12a(3 3 y + Wlfi - 1}, (23) 

where 

Pn = 1 - Pi- (24) 

Now we can write (T231 in the following convenient form: 

. l2aB 3 u(u + u*)(u — u c ) when a^0, 
V"'(M c ,0)={ n ^ P ^ ^ ' T ' (25) 

3/3^yu(yU — /i c ) when a = 0, 
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where 

//* = /i*(a,/3i) = 



l + 6(a,fo) 



(26) 



/x c = /i c (a, ft) = 2r , — Try, (27) 
3#J[1 + 6(a,^i)] 



6 = 6(a - ft) = V 1 + w (28) 

Now as < fit < 1, from (jSJ, we get < 1 - /3 n < 1 ^ > -1 + (3 n > -1 ^ 1 > (3 n > 
< f3 n < 1, consequently, from (128 p . we see that b > 1 for all < /3 n < 1 and for all 
< a < 1 and the sign of equality holds for a = and for all < fi n < 1, i.e., 6=1 only 
when a = and for any admissible value of f3 n . Therefore, from (126]) . we see that /z* is well 
defined as strictly positive real number for all < /3 n < 1 and for all < a < 1 whereas from 
f )2T|) . we see that /z c is well defined as strictly positive real number for all < (3 n < 1 and 
for all < a < 1. Again we note that /x* occurs in the expression of V"'(M C , 0) when a ^ 0, 
i.e., when < a < 1, therefore, in any case (a ^ or a = 0) sign of V'"(M C , 0) depends only 
on the factor \i — fi c , more specifically, V"'(M C , 0) < or V"\M C1 0) > according to fi < \i c 
or > yU c . As /i is restricted by the inequality < < fir < 1 for the case of negatively 
charged dust particles, then the condition /i < fi c <^ V"'(M C , 0) < is always true if /i c > 1. 
But if < n c < n T < 1 then V"'(M C , 0) < or V"'(M C , 0) > according to /j, < /i c or [i > /j, c 
and when /i = /i c , ^/"'(McO) = 0. In this case, i.e., when V'"(M C ,0) = 0, for solitary wave 
solution at M = M c , it is necessary that V""(M c ,0) < 0. So, now our problem, is to find 
the sign of V""(M C , 0) when V"'(M C , 0) = 0, i.e., when /i = /i c . When a^0, with the help 
of simple algebra, we get the following expression of V""(M C , 0) at fi = \i c 

V""(U m 2[(3^-2) 2 (&+l) 2 + 5(& 2 -l)+4&] 

V (Mc ' 0) " 9/33(6 + 1)3 ' ^ 

whereas for a = 0, the expression of V lv (M c , 0) at /i = /i c is given by 



(3/3 n ^2) 2 + l 
9$ 



V" w (M c ,0)= ^ 7 ^ . (30) 



As p n > for all < /3 X < 1, from (ED]), we see that K™(M C , 0) > for a = 0. Again as b > 1 
for a ^ 0, i.e., for all < a < 1 and for all < n < 1, from (ED]), we see that V*"(M C , 0) > 
for a 7^ 0. Therefore, in any case, a = or a ^ 0, for the present problem, we see that for 
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any admissible values of parameters of the system V W (M C , 0) > when /x = /i c , i.e., when 
V"'(M C , 0) = 0. Therefore, there does not exist any solitary wave solution of the energy 
integral <$E§ at M = M c when V"'(M C , 0) = 0. 

Again, differentiating fl25|) with respect to ft, we get 



Now as ft > for all < ft < 1 and p. > 0, from (ED, we see that (d/<9ft)(V w (M c , 0)) < 
for any given value of a, and consequently, V"'(M C , 0) is strictly decreasing function of ft for 
all < ft < 1. Again, it can be easily checked that V"'(M C , 0) > for ft = if /i > /x c (a, 0) 
whereas ^"'(McO) < for ft = if p < p c (a,0). Again, it can be easily checked that 
lim^^i-o V'"(M C , 0) = — /i < 0, where ft — > 1 — <^ ft < 1 and ft approaches to 1, i.e., 
ft approaches to 1 from the left side of 1. Therefore, if /i > /i c (a,0), V"'(M c ,0) strictly 
decreases with increasing ft for all < f3\ < 1 starting from a positive value and ending 
with a negative value whereas if \x < /i c (a,0), V"'(M c ,0) strictly decreases with increasing 
(3i for all < (3% < 1 starting from a negative value and ending with a negative value. 
So, V"'(M C , 0) intersects the axis of (3\ at 0i = (3 C (0 < (3 C < 1) only if /i > /i c (o;,0) and 
consequently, we have the following conclusions regarding the sign of V'"(M c ,0). 

a. If fi > /i c (a,0) and /3 C < f3 1T , we get V"'(M c ,0) > for all < ft < C1 whereas 
V"'(M C , 0) < for all f3 c < ft < ft T and 1/'"(M C , 0) = at ft = ft. 

b. If /j > /i c (a, 0) and ft > ft T , we get V""(M C , 0) > for all < ft < ft T . 

c. If /j < /i c (a, 0), we get V""(M C , 0) < for all < ft < ft T . 

From the above discussions, it is interesting to note that when ft = ft then /i = /i c , since 
for ft = ft, we have ^'"(Mc, 0) = and the equation V"'(M C , 0) = has only one solution 
/i = /i c within the admissible range of fi. On the other hand when /i = /i c then ft = ft, 
since for \i = p c , we have V'"(M C , 0) = and the equation V'"(M C , 0) = has only one 
solution /?i = ft within the admissible range of ft. Therefore, we can conclude that when 
/i = /i c then ft = ft and conversely, when ft = ft then /i = /i c and consequently, at ft = ft, 
we have V""(M c ,0) > 0, i.e., there does not exist any solitary wave solution and/or double 
layer solution of the energy integral at M = M c when ft = ft(^ /i = fi c ). 




when q^0, 



when a = 0. 



(31) 
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So far, we have found that if \i < // c , then only PPSWs are possible at M = M c for all 
< 0i < Pit- On the other hand, if /i > /i c and (3 C < Pit, then there exist NPSW for all 
< fli < (3 C and PPSW for all (3 C < Pi < Pit- However, there do not exist any solitary wave 
or double layer solution at Pi = (3 C . Now, the question is whether Pit is really the upper 
bound of f3\ for the existence of PPSWs. 

For this purpose, consider the existence of PPSW at M = M c for some fii > (3 C and for 
some fixed \i > /i p . Therefore, there exists a (fi m > such that 

V(M C , 0) < for all < <p < </> m , 
V(M c ,(/> m ) = , V'(M c , ( j )m )>0. 

Now as V(M C , <f>) is real for (ft < \1/(M C ) we must have (p m < \I/(M C ), otherwise V(M C , <p m ) is 
not a real quantity. Therefore, 

V{M C , 0) < for all < <p < (<f) m <)^(M C ), (33) 

defines a large amplitude PPSW, which is in conformity with (|32|) . if V(M C , \1/(M C )) = and 

V'(M«0)I*=*(m c ) >0. 

Again let (3i a be the maximum value of Pi up to which solitary wave solution can exist 
at M = M c . As \1/(M) increases with M and M c increases with Pi then we have \1/(M C ) < 
*(M c )\f, 1=f)la . Therefore, 

V{M c ,<f>) < for all < <\> < (0 m < *(M C ) <)*(M c )| A=fto , 

defines the largest amplitude solitary waves if 

^(M c ,*(Af c )| ft=fta = and V'(M C , 0)|^ (Mc)A=/3la > 0. 

Therefore, for the existence of the PPSWs, Pi is restricted by the following inequality: 
Pc < Pi < Pia, where Pi a is the largest positive root of the equation V(M C , ty(M c )\p 1= p la = 
and V'(M C , 4>)\<t>=\s>{M c ),P\=Pi a — f° r au A — ^io- Thus for ji > fi c , PPSW exists in 
Pc < Pi < Pp, where p P = min {/3 la , /3 1T }. 

V. NUMERICAL DISCUSSIONS 

With the help of the above discussions, we have investigated the entire system numerically 
and we have found the following solitary wave solution at M = M c : (i) There exists a value 
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Hp of jj, such that for < fi < fi p there do not exist any solitary wave or double layer 
at M — M c for any value of 0i in < 0\ < Pit- In fact, for any value of fi lies within 
< fi < fi p , there do not exist any non-zero <fi such that V(M C , <fi) = 0; (ii) For any value of 
fi lies within fi p < fi < fi c , there exist only PPSWs for all < 0\ < 0x a ; (iii) In /i c < fi < /V> 
NPSWs exist for < X < C and PPSWs exist in C < X < la ; (iv) For all fi lies within 
/^r < < Ht, NPSWs exist in both < 0\ < 0\ c and 0i c < f3\ < C , whereas PPSWs exist 
for all C < 0i < 0i a . For any point on the curve 0\ = 0\ c in fi r < f> < Ht, one can always 
find a PPDL solution. These findings can be summarized in FIG. [TJ This figure gives the 
exact scenario of the nature of solutions of the system. 

It looks that the solution space is corresponds to nonthermal distributed electrons. How- 
ever, if we move along 0i — axis, we can find different intervals of fi on the basis of the 
polarity of solitary wave solutions in presence of isothermal electrons. In this case, there 
does not exist any solitary wave or double layer for < fi < fi p . Solitary wave starts to exist 
for fi > fi p . Only PPSWs can be found for all ft lies within the interval fi v < fi < fi ci whereas 
for fi lies within fi c < fi < fix , there exist only NPSWs. We have taken two values of /i, one 
lies in ftp < ft < fi c and the other in fi c < fi < fix and draw V(M C , (p) vs. <p in FIG. [2j FIGS. 
[2]^a) and [2](b) verifies the solution space for isothermal electrons. At fi = fi c , there does 
not exist any solitary wave solution, the reasons of which being described in our theoretical 
section. The amplitude of PPSW decreases with increasing fi in fi p < fi < ft c and ultimately 
the PPSWs demolished at fi = fi c . On the other hand, NPSW arrises in the frame after fi 
crosses /i c and the amplitude of NPSW gradually increases with fi in fi c < fi < fi T . Thus 
one may assume that the point fi = fi c acts as a source for NPSWs and as well as a sink 
for PPSWs. It is to be noted that double layer solution cannot be found for Maxwellian 
electrons. 

Now we are to verify the solution space for nonthermal electrons. For this purpose, we 
have taken a value of fi from each of the subintervals fi p < fi < fi c , fi c < fi < fi r , and 
Hr < f- < f-T and draw V(M C , <ft) against <fi in FIGS. [3]-[5]for different values of 0\. In FIG. 
[31 we have taken two values of 0\ lies within < 0\ < (5\ a and both the curves verify the 
existence of PPSW. FIG. H] shows that for a value of fi lies within fi c < fi < fi r there exist 
NPSWs for all < (3 1 < C , whereas only PPSWs exist for all /?i lies within C < X < la . 
Again, (a) and (b) of FIG. [5] show the existence of NPSWs for all < 1 < C and PPSWs 
for all 0i lies within C < 0\ < la , respectively. FIG. [5](c) shows the existence of NPDL 
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solution at ft = ft c for each /i lies within [i r < n < ht- 

A closer look at the solution space given by FIG. [T] reveals an interesting feature near the 
curve ft = ft c . We have seen that NPSWs exist for both ft < ft c and ft > ft c (at least 
in a neighborhood of ft = ft c ) and a NPDL exists at ft = ft c . The immediate question 
is whether the characteristic of NPSWs in < ft < ft c and ft c < ft < ft are same. In 
search of the answer, we draw FIG. [61 In FIGS. [6](a) and[6^b), V(M C , <fi) is plotted against 
<p for three different values of ft, viz., ft c , lc - 0.0005 and ft c + 0.0005. FIG. Mb) shows 
that at ft = ft c , there exists a NPDL whereas, for ft = ft c + 0.0005, there is a negative 
potential soliton. However the curve corresponding to ft = ft c — 0.0005 has no real root 
of (ft in the neighborhood of <fi = <pdi, where <fidi is the amplitude of the NPDL occurs at 
0i — Pic- We notice from figure (a) that V(M C , <ft) again vanishes far beyond <fi = <ftdi for 
ft = ftc, Pi = Pic + 0.0005, and ft = ft c — 0.0005. However, the curves corresponding to 
01 = Pic and ft = ft c + 0.0005 are unable to produce any solitary wave solution but the 
curve corresponding to ft = ft c — 0.0005 gives a NPSW. The amplitude of this solitary wave 
is much higher than that of ft c + 0.0005. Thus there is a jump in amplitude of solitary waves 
and this phenomena has also been observed in some recent works.— ^ Therefore, we can say 
that there exists two types of NPSW, where the first type solitary wave is restricted by 

< ft < ft c and the second type solitary waves exist for all ft c < ft < ft. Moreover, there 
is a jump in amplitude between two types of solitary waves corresponding to the nonthermal 
parameter ft = ft c — e and ft = ft c + e, where e is a sufficiently small positive quantity. 
In FIG. [7] variation in amplitude of NPSWs have been shown for values of ft lies within 
0ic < 0i < 0c whenever fi r < fx < fix- This figure shows that the potential drop is maximum 
at 0i = 0i c and this is reason behind the occurrence of NPDL solution at ft = ft c . This 
figure also shows that amplitude of NPSW decreases with increasing ft and ultimately, these 
NPSWs demolished at ft = ft. FIG. M verifies the fact that there exist NPDL solution at 

01 = 0ic for values of /i in /i r < \i < \it- This figure also shows that the amplitude of 
NPDLs at M = M c decreases with increasing fi in fi r < fi < fi?- Again FIG. [U shows that 
ft c is an increasing function of fi and hence the amplitude of NPDLs at M = M c decreases 
with increasing ft as well. 

In FIG.[9]the variation in amplitude of NPSW and PPSW have been shown for values of ft 
lies within < ft < ft and ft < ft < ft a for jj, c < jj, < ji r . This figure shows that amplitude 
of NPSW decreases with increasing ft and ultimately, these NPSWs demolished at ft = ft, 
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whereas, PPSWs start to occur beyond 0\ = C and amplitude of PPSW increases with 0i 
in C < 0i < 0i a having maximum amplitude at 0\ = 0\ a . Therefore, we can conclude that 
0i = C acts like a sink for NPSWs and source for PPSWs. In fact the same phenomena 
occurs in fi r < \i < /i T when 0i lies within lc < 0i < C an d 0c < 0i < fiia- 

The entire numerical investigation depends on the solution space or compositional para- 
metric space given by FIG. [JJ So for compactness of our investigation, it is customary to 
show that the solution space given by FIG. [1] is unique. To show that this solution space 
is unique we have to show that the values /i c , and fi r maintain the same orderings for 
any value of a as they are in FIG. [TJ To do this, we have numerically investigated these 
values of fi and obtain FIG. [TQJ This figure shows that for any value of a, we always have 
[i p < fi c < f^r- Therefore, for any value of a, the nature of solutions follow FIG. [TJ 

VI. CONCLUSIONS 

The existence of dust ion acoustic solitary wave and double layer at the lowest possible 
Mach number M c have been investigated in an unmagnetized nonthermal plasma consisting 
of negatively charged dust grains, adiabatic positive ions and nonthermal electrons. An 
analytical approach has been constructed to find the polarity of the solitary wave or double 
layer. In the analytical approach a function V'"(M C , = 0) has been constructed, the real 
root 0c of which is a critical value of the nonthermal parameter 0\. Solitary waves exist for 
0i < 0c as well as 0\ > C and the polarity of the solitons have been determined by the 
sign of V"'(M C , = 0). This critical value C of 0\ acts as sink of negative solitons and the 
source of positive solitons. However, it is impossible to find any solitary wave at 0\ = C . 
We have found the upper bounds of the nonthermal parameter for the existence of negative 
and positive potential solitary waves. We have also found the existence of negative potential 
double layer at some point of the nonthermal parameter and observe a jump in amplitude 
of negative solitons at that point, the jump being the amplitude of the double layer formed 
thereat. Depending on these bounds or critical values of the nonthermal parameter, we 
have constructed a compositional parameter space for the system. In this solution space, 
we have delimit the admissible range of fi into three disjoint sub intervals \i v < fi < /i c , 
He < H < fi r , and /i r < \i < fix, by means of three critical points /i c , and /i r of /x. 
However, there do not exists any solution in < /i < fi p . From this solution space we 
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have also found the delimitation of jj, for isothermal electrons. In this case, positive solitons 
exist for ji p < /i < /i c , whereas only positive solitons are possible in /i c < \i < \it- It is 
impossible to find any solitons at /i — /i c . We have not found any double layer solution 
for isothermal distributed electrons. Here, we have considered \xt as the upper bound of 
/i, which is strictly less than 1, since for // — > 1, the effect of charged dust on DIA wave is 
negligible. The amplitude of negative potential solitary waves decreases with increasing (3\ 
while the amplitude of positive potential solitary waves increases with j3±. 
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a = 0.9, M = M 

c 




H -4 



FIG. 1. A graphical presentation of different solitary structures have been given with respect to 
different subintervals of fi within the admissible interval of the nonthermal parameter (3i. In this 
solution space 'P' stands for the existence of PPSW and 'N' stands for the existence of NPSW. 
For any given value of a, NPSWs exist in < fi\ < (3 C and PPSWs exist in j3 c < j3\ < j3\ a - At any 
point on the curve /?i = f3i c , one can always find a NPDL at M = M c . This solution space has 
been drawn for a = 0.9. We have found fi p ps 0.14, /j, c 0.2, and fi r « 0.44. 
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FIG. 2. V(M C ,4>) is plotted against 4> f° r isothermal distribution of electrons with values of ji lies 
within the intervals [i v < fi < fi c (figure (a)) and fi c < /i < (figure (b)). For a = 0.9, we have 
found fip « 0.145 and fi c 0.195. 




FIG. 3. V(M c ,(j)) is plotted against (f> for nonthermal distribution of electrons for values of f3\ lies 
within < j3\ < (i\ a with values of /x lies within the intervals n P < fi < fi c . 
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FIG. 4. V(M C , (f)) is plotted against <fi for two different values of f3\ lies within < f3\ < f3 c (figure 
(a)) and f3 c < /3i < (3\t (figure (b)) with values of fi lies within the intervals fi c < \i < /j, r . For 
a = 0.9, we have found \i c « 0.195, /U r « 0.432 and for = 0.3, we have found /3 C 0.223. 



(a) x10" 5 ( b ) 




«- 4> 

FIG. 5. In (a) and (b), T^(M C ,<^>) is plotted against 4> f° r values of f3\ lies within < f3\ < (3 C and 
/5 C < /?i < /3it> respectively with values of \x = 0.5 lies within the intervals fj, c < /j, < [i r . In (c), 
the same has been plotted at (3\ = (3\ c . For a = 0.9 and for ^ = 0.5, we have found f3 c « 0.427 and 
(3 la « 0.511. 
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FIG. 6. V(M C , 4>) is plotted against 4> for three different values of Pi, viz., Pi c [ — ], (3i c — 0.0005[ ] 

and /?i c + 0.0005[ ] in (a) and (b). In (b) we have shown the region of 4> from (f) = — 2 to <j> = 0, 

whereas in (a) a region of <p from (f> = —11.57 to (j) = —11.54 has been shown. The NPDL profile 
corresponding to f3\ = (3\ c has been shown in (c), whereas the profiles of NPSWs corresponding to 
Pi = Pic — 0.0005 and Pi = Pi c + 0.0005 have been drawn in (d). Profiles in (d) shown a jump in 
amplitude of solitary wave by going from Pi = Pi c — 0.0005 to Pi = Pi c + 0.0005. In all these four 
figures we have used a = 0.9, /(/ = 0.5. 
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a = 0.9, n = 0.5 
1.5 | 




FIG. 7. Variation in amplitude of NPSWs have been shown for values of ft lies within ft c < 
ft < ft. This figure shows that amplitude of NPSW decreases with increasing ft and ultimately, 
these NPSWs demolished at ft = ft. For a = 0.9 and fi = 0.5, we have found f}\ c ps 0.369 and 
P c » 0.427. 




FIG. 8. V(M C , (ft) is plotted against (ft and /?i = j3\ c for four different values of fi, viz., = 0.44 (- 

■ -), n = 0.45 (•••), fJ- = 0.5 ( ) and fi = 0.6 ( ) and respective /?i c 's lies in the neighborhood 

of fa = 0.31001, ft = 0.32203, ft = 0.36917 and ft = 0.43519. 
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FIG. 9. Variation in amplitude of NPSW and PPSW have been shown for values of Pi lies within 
< Pi < [3 C and P c < Pi < Pia- This figure shows that amplitude of NPSW decreases with 
increasing Pi and ultimately, these NPS Ws demolished at Pi = p c , whereas, PPSWs start to occur 
beyond Pi = p c and amplitude of PPSW increases with Pi in p c < Pi < Pi a having maximum 
amplitude at Pi = Pi a . For a = 0.9 and /j, = 0.3, we have found p c « 0.223 and Pi a « 0.33. 
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FIG. 10. fi p , \i r and fi c are plotted against a. From this figure we see that for any value of a, we 
have [ip < fi c < fj, r . 
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